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Abstract 

We study the role that the future detection of the neutrino burst from a 
galactic supernova can play in the reconstruction of the neutrino mass spec- 
trum. We consider all possible Siy mass and flavor spectra which describe the 
solar and atmospheric neutrino data. For each of these spectra we find the 
observable effects of the supernova neutrino conversions both in the matter of 
the star and the earth. We show that studies of the electron neutrino and 
anineutrino spectra as well as observations of the neutral current effects from 
supernova will allow us (i) to identify the solar neutrino solution, (ii) to de- 
termine the type of mass hierarchy (normal or inverted) and (iii) to probe the 
mixing |C/e3p to values as low as 10~^ — 10^'^. 
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1 Introduction 



The reconstruction of the neutrino mass and flavor spectrum is one of the fundamen- 
tal problems of particle physics. It also has important implications for cosmology 
and astrophysics. The knowledge of neutrino masses and mixing will allow us to 
clarify the role of neutrinos in the mechanism of the star explosion and supernova 
nucleosynthesis. 

With the present data on the atmospheric and solar neutrinos, we are taking the 
first steps in the reconstruction of the spectrum. The SuperKamiokande (SK) results 
[H] on atmospheric neutrinos, confirmed by the recent SOUDAN and MACRO 
data, allow us to claim with a high confidence level that the atmospheric neutrinos 
oscillate. Moreover, the oscillations are due to neutrino masses and the mixing in 
vacuum. The data also indicates z/^ ^ z/^ as the dominant mode. All the existing 
experimental results can be well described in terms of the z/^ <-> I'r vacuum oscillations 
with the mass squared difference and the mixing parameters given by ^ 

\Aml^J = (1 - 8) • 10"^ eV^ , sin^ 29 = 0.8 - 1.0 . (1) 

There is no compelling evidence that the electron neutrinos participate in the oscilla- 
tions of atmospheric neutrinos. Moreover, the CHOOZ experiment 0] gives an upper 
bound on the mixing of with Am^ ~ Am^^^: 

sin2 2^e<0.1 for jAm^l > 2 ■ 10"^ eV^ . (2) 

The oscillation interpretation of the atmospheric neutrino data indicates that the 
solution of the solar neutrino problem is also related to nonzero neutrino masses and 
mixing. At the moment, however, there are several possible solutions. Moreover, 
various sorts of data - spectral distortions, day-night effects, seasonal variations - 
favor different possible solutions. A good desription of all the existing data can be 
obtained by 

1. the small mixing angle (SMA) MSW solution: 

Ami = (4 - 10) • 10"^ eV^ , sin^ 29q = (2 - 10) ■ 10"^ , (3) 

2. the large mixing angle (LMA) MSW solution: 

sm^2eQ = 0.7-0.95 , (4) 



sin2 2^0 = 0.8-1.0 . (5) 

e.g. the LOW MSW solution with 
1.0 (see 0, 1^). Results from future 



Ami = (1 - 10) ■ 10"^ ' 
3. the vacuum oscillation (VO) solution: 

2 _ (4 - 6) ■ 10-10 eV^ 
^^0- (6 - 8) • 10-" eV' ' 

Some other possibilities are also not excluded - 
Am^ ~ (0.5 - 2) ■ 10-^ eV2 and sin^ 2^0 = 0.9 - 
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experiments with the existing and new detectors will remove this ambiguity, thus 
identifying the correct solution to the solar neutrino problem. 

Another evidence for the neutrino oscillations follows from the LSND results 0, 
which are not confirmed, but also not excluded by the KARMEN experiment p| . The 
LSND results cannot be reconciled with the solutions of the atmospheric and solar 
neutrino problems in the context of only three known neutrinos, thus requiring the 
introduction of sterile neutrinos. In this paper, we shall consider only the mixing 
between the known three neutrinos (spectra with sterile neutrinos will be discussed 
elsewhere). 

The atmospheric and solar neutrino results (01-0, 0) as well as the existing 
bounds from the other oscillation experiments and the /3/3o,y searches lead to several 
possible spectra of neutrino masses and mixing. The ambiguity is related to (i) the 
unidentified solution of the solar neutrino problem (ii) unknown mixing of z/g in the 
third mass eigenstate which is described by the matrix element Ues (iii) type of hier- 
archy (normal or inverted) which is related to the mass of the third mass eigenstste 
(whether it is the lightest or the heaviest one). The absolute scale of mass is also 
unknown, however this cannot be established from the oscillation phenomena. 



In this paper, we reconsider the effects of oscillations on supernova neutrinos. With 
the existing data on neutrino masses and mixing, we can sharpen the predictions of the 
oscillation effects in the supernova neutrinos. On the other hand, we clarify the extent 
to which the studies of supernova neutrinos can contribute to the reconstruction of the 
neutrino masses and flavor spectrum. We will show that three ambiguities mentioned 
above can in principle be resolved by supernova data. 

The effects of neutrino mixing on the neutrino fluxes from the supernova have 
been extensively discussed in the context of 2z/ mixing. For a wide range of mixing 
parameters (Am^ < 10^ eV^), the neutrinos encounter their MSW resonance densities 
inside the star, hence the studies of resonant neutrino conversions inside the star [p!0[ - 
1^ are crucial. For very low values of Am^ (e.f/. Am^ < 10~^^ eV^), the vacuum 



oscillations on the way from the star to the earth need to be taken into account |T^, T^\. 
In the presence of a strong magnetic field, the spin- flip effects become important ||15|| : 
the spin-flavor precession and resonant spin-flavor conversions |jT^ may affect the 
observed neutrino fluxes. If sterile neutrinos are involved in the neutrino conversions, 



they may enable the r-process nucleosynthesis []18 



The effects of the neutrino conversions can be observed through, e.g. (i) the 
disappearance (partial or complete) of the neutronization peak; (ii) the interchange 
of original spectra and the appearance of a hard z/g spectrum; (iii) distortions of the 
Pf. energy spectrum; (iv) modification of the spectrum (in particular, the effects of 
large lepton mixing on the Ue spectrum have been extensively studied |1^); (v) earth 
matter effects. The observation of the neutrino burst from SN1987A [20| has already 
given bounds on the large mixing of active neutrinos (|TT|, |T^,|T^-||2^) and on the 



mixing of z/g with sterile neutrinos |23 



The main features of transitions of supernova neutrinos in the case of mixing 
(p2[,[^-[0]) can be understood in terms of the 2z/ mixing. The system has two 
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resonances^ Under the assumptions of mass hierarchy and smallness of mixing, the 
dynamics of the two level crossings splits. As a result, the factorization of probabili- 

m 



ties occurs 
place 



25 



In the presence of sterile neutrinos, multi-level conversions take 
which may be interpreted in terms of the constituent 2z/ conversions 
for small mixing angles and Am^ hierarchy M. 
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In this paper, we study the conversions of supernova neutrinos in the Su context, 
taking into account the recent results on the neutrino masses and mixing. We consider 
the effects for all possible schemes of neutrino masses and mixing which explain the 
atmospheric and solar neutrino data. For each scheme, we find the modifications of (i) 
the neutronization peak (ii) the energy spectrum and (iii) the energy spectrum, 
which can be observed directly. We also determine the spectrum of the non-electron 
neutrinos which can in principle be studied by neutral current interactions in reactions 
with different energy thresholds. 

The paper is organized as follows. In Sec. II we describe the features of initial 
neutrino fluxes from the supernova and the dynamics of neutrino conversion on their 
way out to the surface of the star. In Sec. Ill, we derive general expressions for the 
transition probabilities for the schemes with normal mass hierarchy. We also calculate 
the earth matter effects on the neutrino spectra. In Sec. IV, we find the final neutrino 
spectra at detectors for the schemes with the normal mass hierarchy. In Sec. V, we 
perform similar studies for the schemes with the inverted mass hierarchy. In Sec. 
VI, we discuss the observable signals and the signatures of various mixing schemes 
In Sec. VII, comparing results for various schemes we conclude about the possibility 
to discriminate the schemes by future observations of neutrino bursts from a galactic 
supernova. 



2 Mass spectra, fluxes and dynamics of conversion 

In this section, the generic properties of the initial neutrino fluxes will be summarized. 
We identify the neutrino mass and mixing parameters relevant for the supernova neu- 
trino conversions, and consider main aspects of dynamics of neutrino conversion inside 
the star: the transition regions, the factorization of dynamics and adiabaticity. Fi- 
nally, we construct the level crossing schemes for the normal and the inverted mass 
hierarchy. 



2.1 Neutrino fluxes 

In what follows we will summarize the generic features of the original fluxes which do 
not depend on the model and the parameters of the star. The deviations from these 
features will testify for neutrino conversions. 

•^The radiative corrections to m^^ and m^^ imply the existence of one more resonance between 
the two non-electron neutrinos, but since the two non-electron neutrinos cannot be distinguished at 
the detector, the conversions between them do not affect the observations. See sec. |2.6|. 
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1. Flavor of the neutronization peak. During first few milliseconds of the 
neutrino burst from a supernova, the signal is expected to be dominated by the z/g, 
which are produced by the electron capture on protons and nuclei while the shock 



wave passes through the neutrinosphere [^. Since the original flux is z/g, the final 



observed fluxes give a direct measurement of the extent of conversion of v^. into the 
other neutrino species. 

2. Inequalities of average energies of the spectra. Since v^, interact more 
strongly with matter than the other species, their effective neutrinosphere is outside 
the neutrinospheres of the other species and hence they have a lower average energy 
than i/e and v^. The v^. also interact via charged current, but the cross section is 
smaller, so their average energy is more than that of the z/e, but less than that of v^- 

«)<«)<«) • (6) 

Here subscript refers to the original spectra. The neutrino conversion changes in- 
equalities (|]). 

3. The pinched spectra during the cooling stage. Let and Fg be the 

original fluxes of z/g and Ve respectively, produced at the cooling stage. The "non- 
electron" neutrinos (z/^, z/^, z/^, z>r) have the same neutral current interactions inside 
the supernova, and their original fluxes are expected to be approximately equalQ. 
In what follows we will neglect the difference of fluxes. We will denote these four 
neutrino species collectively as and the original flux of each of them as 

The spectra of neutrinos from the cooling stage are not exactly thermal. Since 
the neutrino interaction cross-section increases with energy, even for the same species 
of neutrinos, the effective neutrinosphere radius increases with increasing energies. 
Then, even if the neutrinos at their respective neutrinospheres are in thermal equi- 
librium with the matter, the spectrum gets "pinched" , i.e. depleted at the higher as 
well as the lower end of energies in comparison with thermal spectrum. (See Fig. |^ 
for an illustration.) It can be shown that pinching of the instantaneous spectrum is 
the consequence of the following two facts: (i) the temperature inside the supernova 



decreases with increasing radius and (ii) the density decreases faster than 1/r p4 |. 
The pinching of the instantaneous spectrum can be extended to the pinching of the 
time-integrated spectrum as long as the time variation of the average energy of the 
spectrum is small. Numerical simulations of the neutrino spectra confirm the pinching 
of spectra [^, 

One way to parametrize the pinched neutrino spectra is to introduce an effective 
temperature Tj and an effective degeneracy parameter 77^ (which has the same sign 
for neutrinos and antineutrinos and cannot be considered as the chemical potential) 

^ The presence of real muons in the central part of the star leads to a nonzero chemical potential 
of the muon neutrinos and hence to a difference of fluxes |^3| . However, in the neutrinosphere with 
T — % MeV, the concentration of muons is smaller than 1%. 
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Figure 1: The number of u — N charged current events for (a) a thermal spectrum 
(T = 3, ?7 = 0) (b) a "pinched" spectrum (T = 3, r/ = 3). T and r] are the parameters 
as in Eq. (|^). Both the energy resolution function and the lower energy threshold 
effects have been taken into account. The spectra are normalized to have equal areas. 



in the Fermi-Dirac thermal spectrum for each species i: 

For a pinched spectrum, rji > 0. The value of rji is the same for all species (neutrinos 
as well as antineutrinos, since they have the same interactions), and are in general 
different from rje or r]g. The value of 7]i need not be constant throughout the cooling 
stage. Typically, 

^ 3 - 4 MeV , ^ 5 - 6 MeV , ^ 7 - 9 MeV . (8) 



The typical values of rji are ^ 



ri,^3-5 , r/e-^2.0-2.5 , r/, ^ - 2 . (9) 

Notice that the strongest pinching is expected for the z/g spectrum. These values, 
however, are model-dependent and we shall use them only as a guide. 

In what follows, we shall calculate the fluxes of electron neutrinos F^, electron an- 
tineutrinos, Fg, and the total flux of the non-electron neutrinos, "4Fa;" at the Earth 
detectors. 



2.2 Neutrino mass spectra 

We consider the system of three active neutrinos Uf = (z/e, f^, Ur) mixed in vacuum, 
such that 

Vf = Uv , (10) 
6 



where z7 = (z/i, z/2, z/3) is the vector of mass eigenstates and 

U^WUf^W , / = e,/i,r, , ^ = 1,2,3 (11) 

is the mixing matrix. We take 

Am^2 = ml-ml = Aml^^ , (12) 

so that the oscillations driven by Am'^^ solve the atmospheric neutrino problem. We 
identify 

Amji = ml — ml = Am^ , (13) 

where Am^ is in one of the regions (SMA, LMA or VO) implied by the solar neutrino 
data (eqs. |,|,D. 

The key features of the spectra which play an important role in the applications 
to the supernova neutrinos are: 

1. the hierarchy of Am^: 

|Am^2| ^ lAm^^l > \Aml^\ ; (14) 

2. the upper bound on Am^: 

\Aml\ < 10-2 eV^ . (15) 

Since and Ur are indistinguishable in the supernova neutrino studies, the neu- 
trino transitions are determined by the mixings of the electron neutrino only, i.e. by 
the elements Uei (0) (see sect. ^ for more details). Moreover, the three Uei^s are 
related by the unitarity condition J2 l^ejP = 1, so that only two mixing elements are 
relevant, and one can use |[/e2| and \Ue3\. 

The element t/es is small, as mentioned in the introduction. If |t/e3p 1, then 
Ue2 (and therefore, t/gi) can be found from the solar neutrino data: 

MUe2\^\Uel\^-MUe2\\l-\Ue2\^)-Sm^2eQ , (16) 

where Oq is the mixing angle determined in the 2i> analysis of the solar neutrino data. 

The system is then determined by two pairs of parameters (Amf , sin^ 26'j), i = 
L,H, where 

(Ami,sin2 2^L) ^ (Am|, sin^ 2^0) 

{Aml,sm'2eH) ^ {Aml^AlUesl') • (17) 

Correspondingly, it can then be described by two points in the (Am^, sin^ 2^)-plane. 

The current oscillation data do not determine the mass and flavor spectrum com- 
pletely. As we have already mentioned in the introduction, the uncertainty is related 
to: 

1. The discrete ambiguity in the solution of the solar neutrino problem: The data 
favor three solutions indicated in the introduction (SMA, LMA and VO), and 
some other solutions {e.g. the low Am^ MSW solution ||^, ^, the trimaximal 
mixing solution |Q), although disfavored, are not excluded. The future solar 
neutrino experiments will remove this ambiguity and sharpen the determination 
of the oscillation parameters. 
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2. The ambiguity in the sign of Am^g (and Am|^): this determines the type of 
neutrino mass hierarchy. We refer to the case 

Am32 > , i.e. m3>m2,mi (18) 
as to the spectrum with normal mass hierarchy and to the case 

Am32 < , i.e. m2,mi>m3 (19) 
as to the spectrum with inverted mass hierarchy. 

The key difference between these two hierarchies is that, in the normal hierarchy, 
the small t/gs-admixture of is in the heaviest state whereas in the inverted 
hierarchy, this admixture is in the lightest state. 

The type of hierarchy can in principle be established in future studies of atmo- 
spheric neutrinos. If the effects of u^. oscillations will be observed and the sign 
of the charged lepton produced by the atmospheric neutrinos will be identified 
in the future experiments, the studies of matter effects in the neutrino and an- 
tineutrino channels will allow to establish the sign of Am^. The sign can also be 
found from the studies of matter effects in the future long baseline experiments, 
in particular, with neutrino factories. 

3. The value of f/es is unknown. In principle, future atmospheric neutrino exper- 
iments and the long baseline experiments will be able to measure or further 
restrict Ues- 

4. In the case of the SMA and LMA solutions, Am^i > 0. But the sign of Arri^i 
is undetermined in the VO solution. The feasibility of resolving this ambiguity 



has been recently discussed in ||39 | 



Summarizing, the ambiguity in the present analysis is related to the solution of 
the solar neutrino problem, the type of hierarchy and the value of Ues- The first two 
ambiguities lead to six possible schemes of neutrino masses and mixing. Within each 
scheme, the predictions depend on the value of f/es- 

2.3 Neutrino conversion regions 

In supernova, the transitions occur mainly in the resonance layers, where the density 
varies between (p^es — ^Pres) (Pres + ^Pres)- Here pres is the resonance density: 

^ cos20 , (20) 



^ 2^Gp E Y, 

Gp is the Fermi constant, m^r is the mass of the nucleon, E is the neutrino energy 
and is the electron fraction - the number of electrons per nucleon^. For small 

^It is worthwhile to introduce the resonance density for "VO" parameters too, ahhough the 
dominating effect could be the vacuum oscillations. Wc shall consider the limit of sin^ 2(? ^ 1 as a 
special case. 
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vacuum mixing, the width of the resonance layer equals 

2/\pres ~ 2p^e^tan26' . (21) 

Using (POD, the resonance matter density can be written as 

, / Am^ \ /lO MeV\ /0.5\ ^ , , 

Pres - 1.4 X 10*^ g/cc [y^] [-^—) ( ^°«2e ■ (22) 

There are two resonance layers. The layer at higher densities (if-resonance layer), 
which corresponds to Am^^^, is at 

~ 10=^ - 10' g/cc , (23) 

and the layer at lower densities (L-resonance layer), characterised by AttIq, is at 

r 5-15 g/cc for SMA 
Pl=\ 10-30 g/cc for LMA (24) 
[ < 10-4 g/cc for VO . 

The regions where the neutrino transitions occur are thus far outside the core of 
the star - in the outer layers of the mantle. Therefore 

• the transitions do not influence the dynamics of collapse or the cooling of the 

core; 

• the r-processes, which occur at p > 10^ — 10^ g/cc, are also not affected; 

• the shock wave does not influence the neutrino conversions (indeed, during the 
time of cooling by neutrino emission {t ~ 10 sec), the shock wave can only reach 
layers with densities p > 10^ g/cc [^, p[]); 

• the density profile encountered by the neutrinos during their resonant conver- 
sions is almost static, and the same as that of the progenitor star. 

In the region with densities > 1 g/cc, the electron fraction is almost constant and 
the density profile can be approximated by [|1^, ^ 



pre~2-104g/cc — , forp>lg/cc. (25) 

VIO^ cm / 

For p < 1 g/cc, the fraction of hydrogen increases and Ye becomes larger than 0.5. 
The exact shape of the density profile depends on the details of the composition of 
the star. 
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2.4 Factorization of dynamics 

The hierarchy of Am^, and therefore the hierarchy of the densities of the resonance 
layers, leads to the "factorization" of dynamics of conversion: the transitions in the 
two resonance layers can be considered independently and each transition is reduced to 
a two neutrino problem. Indeed, in the if-resonance region, the mixing associated 
with AttIq is suppressed by matter. The suppression factor is 

^e2 PL<lQ-2 . (26) 



r\-/ 



Ue2 PH 

Correspondingly, the effects driven by AitLq are suppressed by more than two orders 
of magnitude. 

In the L-resonance region, the mixing associated with Am^^^ coincides with the 
vacuum mixing: ~ Ues- The matter corrections are strongly suppressed: 

C/- = [/e3(l + O(0) , — <10-' . (27) 

Ph 

That is, the mixing associated with Am^^^ is almost constant, and therefore the level 
z/3 practically does not participate in the dynamics. By an appropriate redefinition of 
the fields, the problem can be reduced to a two state problem. The state 1/3 decouples 
from the rest of the system, producing just an averaged oscillation effect [pG] . 



If the mixing Ue^ is very small, the resulting survival probability of i/g (^e) is also 



factorized 24, 25 : 



[Pil^e^l^e)]total = [Pil^e^l^e)]H X [P(zy, ^ zyj]^ +0(|f/e3|') , (28) 

and similarly for u^- 

2.5 Adiabatic, non-adiabatic and transition regions 

The dynamics of transitions in each resonance layer is determined by the adiabaticity 



parameter 7 |]43|, ^ 



_ Am^ sin^ 2^ 1 
'^^ 2E cos 2^ {l/ne){dnjdr) ' ^^^^ 

such that the "flip" probability - the probability that a neutrino in one matter eigen- 
state jumps to the other matter eigenstate - is 

Pf = exp (30) 



as given by the Landau-Zener formula |^3|. Adiabatic conversion corresponds to 
7^1, i.e. to a very small flip probability Q. 

^The Landau-Zener formula is valid for a linear variation of density in the resonance region and 
a small mixing angle 9. For an arbitrary density distribution and mixing angle, the Landau-Zener 
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Let us consider the density profile of tlie form 

p = A/r" , (32) 



where A is the proportionahty constant. From (|29|) , we get the adiabaticity parameter 
for this profile as 

^ J_ (Anf\ sin^ 26 (2^20^ \ 

^ 2n[ E J (cos 2^)1+1/" [ mj^ J " ^ ^ 

Here we have used the resonance condition to express r through the oscillation pa- 
rameters. Note that the dependence of 7 on the absolute scale of density, A, is rather 
weak: 

For n = 3, the change of A by one order of magnitude leads to the change of 7 by a 
factor of 2. 

For a fixed density scale A, the value of 7 depends on the power index n as 

7 oc 1/n , (35) 

so that a variation of n between 2 to 4 leads to the variation in 7 by a factor of 2. 

In turn, the uncertainty of a factor of 2 in 7 is equivalent to the change of sin^ 29 by 
a factor of 2 (at small values of 9), or the change of Am^ by a factor of 2"/*^""^^ ~ 3—4. 

The lines of equal 7 (and therefore, equal Pj) on the (Am^ — sin^ 29) - plot are 
determined by 

(Am^)^-^/" sin^ 29 = const (36) 

for small 9. 

For n = 3 which will be used in the further calculations, we get 

Am^ = const / sin^ 29 . (37) 
From (|30|) and (|33D , the flip probability as a function of energy can be written as |2^ 

-E 

-'-'na 1 



Pf = Exp 
where 



E 



(St 



1/2 



^„„.(M-^!I^!|^(!^,, . ,39) 

\12/ cos^2t^ Y m^v / 

The dependence of Pf on E/ E^a is shown in Fig. |^. One can divide the whole range 



formula (pQ) gets modified to 



exp{~^^F) ~ ea;p(-f 7F/sin^ 



where is a function of the density profile and the mixing angle. The adiabaticity condition is, 
approximately, 7^" ^ 1. At small mixing angles and linear density variation in the resonance region, 
F « 1 and (pll) reduces to (|^. 
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of energy in three parts: 

For E/Ena < 10^^ (region I), we get Pf ~ 0. In this range, pure adiabatic 
conversion occurs. 

For E/Ena > 10^ (region III), the flip probability is close to 1, which corresponds 
to a strong violation of adiabaticity. 

In the transition region E/Ena = 10^^ — 10^ (region II), increases with the 
neutrino energy. This region spans almost three orders of magnitude in energy, which 
is substantially larger than the range of energies in the neutrino spectrum. Notice 
that for the exponential density distribution (which is the case inside the sun), the 
transition region is narrower (about two orders of magnitude) and correspondingly, 
the energy dependence in this region is stronger. 

The observable part of the supernova neutrino spectrum lies mainly between the 
energies of 5 and 50 MeV, i.e. it spans about one order of magnitude. If the spectrum 
is in region I, completely adiabatic conversion occurs for the whole spectrum. In the 
region II, the conversion depends on energy, however the dependence is not strong 
over the relevant range of energies. The average energies of any two neutrino species 
differ by a factor of less than 3, and for AE / E ~ 3, the variation in the flip probability 
is APf < 0.2 (from Fig. In the first approximation, the final spectrum can then be 
characterized by an average or effective flip probability. This is illustrated in Fig. |^: 
the spectra of the number of events taking into account the energy dependence of Pf 
are shown in (a) and the spectra with an effective flip probability (Pf) are shown in 
(b). It may be observed that the spectra with an appropriate value of (Pf) can mimic 
the features of the actual spectra. 

In Fig. 1^, we show the contours of equal flip probability Pf (|3T|) in the (Am^ — 
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Figure 3: The number of Uf, — N charged current events, taking into account the 
energy dependence of the flip probabihty. The parameters for the original spectra are 
taken to be Tg = 3 MeV, ?7e = 3, = 8 MeV, rj^ = 1. In (a), the flip probability is 
Pf = Exp[-{E^a/E)^/^] with (al): = 0.05, (a2): E^a = 2, (a3): E^a = 10, (a4): 
Ena = 50 MeV. In (b), the effective flip probability (Pf) is (bl): 1.0, (b2): 0.85, (b3): 
0.6, (b4): 0.0 



sin^ 29) plot for two different energies on the borders of the observable spectrum. We 
also show the parameter ranges which explain the solar and atmospheric neutrino 
data. The dark band (the "atmospheric neutrino band") corresponds to the allowed 
range of Amgj^. The rightmost part of this band is excluded by the CHOOZ exper- 
iment 0]. The range of |[/e3p that can be probed by the long baseline experiment 
MINOS [|§ is also shown. 

The contours of Pf = 0.1 and Pf = 0.9 divide the plot into three regions, corre- 
sponding to the three regions in fig. |^: 

• The "adiabatic region" (I) is the region above the contour with Pf = 0.1, where 
adiabaticity is well satisfied and strong flavor conversions occur. 

• The "transition region" (II) is the region between Pf = 0.1 and Pf = 0.9 
contours. Here the adiabaticity is partially broken and the transitions are not 
complete. Moreover, the extent of transitions depends on the energy. 

• The "non-adiabatic" region (III) lies below the Pf = 0.9 contour. The neutrino 
conversions are practically absent. 

Since the adiabaticity breaking increases with E and decreases with the increase of 
Am^ and sin^ 29, the lines of equal Pf for E = 50 MeV are shifted to larger Am^ and 
sin^ 29 relative to the lines for E = 5 MeV. The dependence of the contours of equal 
Pf on the density profile is rather weak, as can be seen from the previous discussion. 
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Figure 4: The contours of equal flip probability Pf. The solid lines denote the 
contours of flip probability for a 5 MeV neutrino: the line on the left stands for 
Pf = 0.9 (highly non-adiabatic transition) and the line on the right stands for Pf = 0.1 
(adiabatic transition). The dashed lines represent the corresponding flip probabilities 
for neutrinos with energy 50 MeV. SMA, LMA and VO correspond to the solutions 
of the solar neutrino anomaly. The two vertical lines indicate the values of 4|?7e3p = 
sin^ 29 lying on the borders of the adiabatic, non-adiabatic and transition regions for 
Am^ corresponding to the best flt value of the atmospheric neutrino solution. 



Even with conservative estimates, the borders of the regions have an uncertainty of 
a factor of 2 in sin^ 29 for a given value of Am^. 

As follows from Fig. ^, the LMA solution lies in the adiabatic region iQ. The SMA 
solution is in the transition region. The VO solutions are either in the transition 
region or in the non-adiabatic region, which depends essentially on the density proflle 
in the outermost layers of the star (p < 1 g/cc) and the precise value of Am^. 

As was described in sect. |2.2| , each neutrino mass and flavor spectrum can be 
represented by two points (|17D in the (Am^,sin^2^) plot (Fig. One point, cor- 
responding to (Amgj^, sin^ 2^e3), should lie in the atmospheric neutrino band, and 
the other point corresponding to (Am|j^, sin^ 2^e2) should lie in one of the "islands" 
corresponding to the solutions of the solar neutrino problem. These two points char- 
acterize the layers H and L (|^,|^) respectively. Let Ph{Ph) and Pl{Pl) be the 
probabilities that the neutrinos (antineutrinos) jump to another matter eigenstate in 
these two layers. The extent of conversion is determined by the values of these four 
flip probabilities. 

From Fig. ^, we conclude that the if -resonance is in the adiabatic range (region 

^The LOW solution also lies in the adiabatic region, so all the results for the LMA scenario are 
also valid for the LOW scenario. 
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I) for 

sin2 2e,3 = 4|[/e3p>10-3 , (40) 
and in the transition region (region II) for 

sin^ 2^e3 ~ 10"^ - 10"^ . (41) 

As we shall see in sect. H and ^, the features of the final spectra depend strongly on 
the region in which the H resonance lies. According to (0), the supernova neutrino 
spectra are sensitive to as low values of 4|?7e3p as 10^^ — 10^^. This is more than two 
orders of magnitude better than the current bounds [Q or those expected from the 
planned long baseline experiments. 



2.6 The level crossing schemes and initial conditions 

In the basis of flavor eigenstates (z/g, z/^, z/^), the evolution of neutrinos at densities p < 
10^ g/cc relevant for neutrino conversion (see below) is described by a Schrodinger-like 
equation with the effective Hamiltonian 

H = + V 

2E 



1 

2E 



( ml^ + 2EV 



(42) 



er 




where V ~ Diag{V, 0,0), and V = y^Gpne is the effective potential for the electron 
neutrinos due to their charged current interactions with electrons. 

Since any rotation in the (z/^ — z/^) subspace does not affect the physics, it is 
convenient to perform a rotation of the neutrino states (z/g, z/^, z^r) ~^ {i^eii^fi,'ii^T')i 
which diagonalizes the (z/^, z/^) submatrix of ( ^2]) |]47[] . (The potential V appears 
only in the element Hee, and hence is not affected by this rotation.) The effective 
Hamiltonian in the new basis becomes 




ml^ + 2EV ml^, ml^, 

^=7^1 <' ^'v I • (43) 

mi 



T T 



At V ^ my/{2E), the off-diagonal terms can be neglected and the Hamiltonian 
3]) becomes diagonal: 

H ^DmgiV, ml, ^„ml,^,) . (44) 

That is, the basis states (z/g, z/^/, z/^/) are the matter eigenstates. These are the states 
that arrive at the conversion regions as independent (incoherent) states and transform 
in this region independently. 

Notice that a difference between the potentials of and z/^ appears in the second 
order in the weak interactions due to difference of masses of the fi and r charged 



leptons EH 



V,. ^ V^£^ In^ - 1 + 10- V , (45) 
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where is the r mass, mw is the VT-boson mass, and Yn are the numbers of 
electrons and the neutrons per nucleon respectively. Therefore a complete form of 
the matrix of potentials is V ~ DiagiV, 0, V^t)- 

The potential Vfj_r becomes important at high densities. One has 

V ~ ^rnlJ2E ^ 2mlJ2E at p,, ~ 10^ - 10« g/cc. 

At p ^ p^ri and in particular in the region of the neutrinosphere, the potentials 
V and V^r dominate over the other terms in the Hamiltonian, and the Hamiltonian 
becomes approximately diagonal: if ^ V ~ Diag(\^, 0, V^r)- This means that at high 
densities the flavor states coincide with the eigenstates in medium. 

Let us recall that the non-electron neutrinos are produced in the neutral current 
processes which are flavor blind, i.e. they produce a coherent mixture of matter 
eigenstates. This coherence, however, disappears in the evolution that follows. The 
Vf. departs from the coherent state due to the large potential V , whereas the coherence 
of and is broken by V^t- 

In the interval of densities where V^r ^ ^'^'itm/'^^ ^ ^^e potential V^^ can be 
neglected, so that we arrive at the Hamiltonian ( ^3]) with eigenstastes {u^, t^^', i't')- In 
the region V^^ ~ Aml^^/2E, the level crossing occurs which leads to the transitions 
z/^ — > u^j. and z/^ — > u'^. In the antineutrino channel, we have z/^ —>■ u'^ and Ur —>■ v'^- 

Since the initial fluxes of and are equal we get that the fluxes of z/^/ 

and Vr' will be also equal Therefore {l'e^^^l' i^-'t') with fluxes {F^,F^,F^) and 

correspondingly (t'e, J^^t', ^^r') with fluxes {Fg,F^,F^) can be considered as the initial 
state in our task. 



The Hamiltonian p3|) allows us easily to construct the level crossing scheme. In 
Fig. H, we show the generic level crossing diagrams for the normal and inverted mass 
hierarchies, for small 6q (SMA) as well as large ^0 (LMA or VO). The diagonal 
elements of the Hamiltonian (^), Hii{ne) {i = e, yu',r'), determine the energies of 
the flavor states shown by the dotted lines. The crossing of these levels indicates a 
resonance. H and L are the two resonances at higher and lower densities respectively 
The solid lines represent the eigenvalues of the Hamiltonian (^). 

In the case of antineutrinos, the effective potential V for the z/g has the opposite 
sign : V = —y/^Gpne. The antineutrinos can then be represented on the same level 
crossing diagram, as neutrinos travelling through matter with "effectively" negative 
rig. The half-plane with positive values of Ue then describes neutrinos and the half- 
plane with negative values of rie describes antineutrinos. 

The neutrinos (antineutrinos) are produced inside a supernova in regions of high 
matter density. On their way towards the earth, they travel through a medium with 
almost monotonically decreasing density, towards the vacuum where both neutrinos 
and antineutrinos have vanishing effective potentials. This corresponds to starting at 
the right (left) extreme ends of the Ue axis in Fig. ^ and moving towards = 0. 

The H resonance lies in the neutrino channel for the normal hierarchy and in the 
antineutrino channel for the inverted hierarchy. The L resonance lies in the neutrino 
channel for both the hierarchies as long as the solar neutrino solution is SMA or LMA. 
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Figure 5: The level crossing diagrams for (a) the normal mass hierarchy and small 
Qq^ (b) the inverted mass hierarchy and small ^0, (c) the normal mass hierarchy 
and large ^^0, (d) the inverted mass hierarchy and large ^0. Solid lines show the 
eigenvalues of the effective Hamiltonian as functions of the electron number density. 
The dashed lines correspond to energies of flavor levels z/g, v^i^ and v^i. The part of 
the plot with < corresponds to the antineutrino channel. 
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For the VO solution, the L resonance may he in either of the two channels, neutrinos 
or antineutrinos. 



3 Conversion probabilities and the neutrino fluxes 
at the detectors. The case of normal mass hier- 
archy 

In this section, we derive general expressions for the transition probabilities using the 
level crossing scheme for the normal mass hierarchy as shown in Figs, ^a andj^c. (The 
inverted hierarchy will be discussed separately in Sec. ||.) 

3.1 Probabilities of conversion inside the star 

As has been discussed in sec. p.6| , the neutrino fluxes arise from the central part of 
the star, in the region of high density. For p ^ PhiPli where all the mixings are 
highly suppressed, the flavor states {ve^^^i' i^-'t') coincide with the eigenstates in the 
medium : 

'^Sm = T^e , ^2m = ^r' , ^Im = T^^i' ■ (46) 

The original fluxes of neutrino eigenstates in the medium equal 

Let us calculate the fluxes of the mass eigenstates z/j at the surface of the star. 
These states, being the eigenstates of the Hamiltonian in vacuum, travel indepen- 
dently to the surface of the earth. 

Taking into account that the dynamics of transitions in the two resonance layers 



are independent (see Sec. ^^) , the fluxes of neutrino mass eigenstates at the surface of 
the star can be written down directly by tracing the path of the neutrinos in the level 
crossing diagram. We find the modifications of fluxes in terms of the flip probabilities 



Ph, Pl, Ph and Pl introduced in Sec. 2.5 



Let us first calculate the flux of vi at the surface of the star. There are three 
independent contributions to this flux, from the initial Ve-, ^ii' and fluxes. The 
probability that the original state Ve (which coincides with v^^n in the production 
region) arrives at the surface of the star as vi is PhPli since the state has to flip to 
the other matter eigenstate at both the resonances. The contribution to the final vi 
flux from the original Ve flux is then PhPlF^. Similarly, the contributions from the 
original i/^/ and z/,-/ equal (1 — Pl)F^ and Pl(1 — Ph)F^ respectively. The total z/i 
flux at the surface of the star equals the sum of the three contributions: 

F^=PhPlF^,+{1-PhPl)F^, . (48) 

Similarly, the fluxes of neutrino mass eigenstates V2 and z/3 arriving at the surface of 
the star are 

F2 = {Ph-PhPl)F! + {1-Ph + PhPl)F^. , 

F3 = il-PH)F^ + PHF^^ . (49) 
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Due to the divergence of the wavepackets, any coherence between the mass eigen- 
states is lost on the way to the earth. Indeed, over a distance L, the two wavepackets 
corresponding to two mass eigenstates with a given Am^ and having an energy E 
separate from each other by a distance 



AL 



2^2 



(50) 



Even for the smallest Am^ ~ 10"^° eV^, for ^ ~ 10 MeV and L ~ 10 kpc ~ 10^^ 
cm, we get AL ~ 10~^ cm. The lengths of the individual wavepackets are a < 1/T ~ 
10~^^ cm (where T is the temperarture of the production region) is much smaller. 

The spread of the wavepackets implies that the neutrinos arrive at the surface of 
the earth as incoherent fluxes of the mass eigenstates. Upto a geometrical factor of 
l/L^, they coincide with the fluxes given in (^,|^). We can rewrite them as 



Fi = ttiF^ + {1 - ai)F^ , 



with 



ai = PhPl, a2 = Ph{1 - Pl), = I - Ph ■ 

The factor of l/L^ is implicit in the fluxes at the earth. 

Using (|5T|), we find the net flux of electron neutrinos at the earth: 

Fe ^ ^ I Uei I Fi 

i 

= F^J2\Ue4\ + F^{l-J2\Ue.\'a,) , 



(51) 



(52) 



(53) 



where we have taken into account the unitarity condition J2 l^eiP 
electron neutrino flux reaching the earth can thus be written as 



Fe=pF° + (l-p)F° 



1. The final 



(54) 



where 



P 



^ ^ I Ugi I di 

i 

\Uel?PHPL 



\Ue2? 



H 



PhPl] 



\U,:,WI-Ph) 



(55) 



According to (|5^), p may be interpreted as the total survival probability of electron 
neutrinos. 

The original total flux of the neutrinos z/g, z/^, Vr is F^+2F^. Using the conservation 
of flux, we find the combined flux of and Vr at the earth (F^ + Fr) as 



F^ + F. = (l-p)F° + (l+p)F° 



(56) 



Note that the final fluxes of the flavor states at the earth ( |5^ , |5BD can be written 
in terms of only the survival probability p. This is a consequence of two facts: 1) 
at each transition, one of the neutrinos is decoupled so that the task reduces to 2v 
mixing, and 2) the original fluxes of v^i and z/,-/ are equal. 
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3.2 Conversion probabilities for antineutrinos 

Let us consider the antineutrino transitions. In the high matter density region 
(p ^ Ph,Pl), the antineutrino flavor eigenstates coincide with the eigenstates in 
the medium as (see Figs. |c): 

i^lm = t'e , i^2m = t'/x' , ^'Srn = J^t' , (57) 

SO that the original fluxes of antineutrino eigenstates in the medium equal 

PO _ pO _ pO _ 77.0 /c-oN 



The small mixing angle 6e3 is further suppressed in the medium, so the z/g ^ z/3 
transitions are negligible. The state usm, being far from the level crossings, propagates 
adiabatically: Ur' 1^3 ■ Depending on the parameters of the solution for the solar 
neutrino problem, the propagation of the other two states may be adiabatic or non- 
adiabatic |T^. From considerations similar to those in the neutrino channel, we get 



f, = pf^ + {i-p)f: 

where p, the survival probability of Ue, equals 





X 5 



P 



\Uel\^(l-PL) + \Ue2\^PL 



In the case of completely adiabatic propagation, = and 

P=\Uel\'. 



(59) 



(60) 



(61) 



Let Ffi + Ff be the combined flux of z/^ and Ur. From (59) and the conservation 
of flux, we get the combined flux of non-electron antineutrinos: 

F^ + F, = (l-p)F°+(l+p)F° . (62) 

the total flux of the non-electron neutrinos (including an- 



(63) 



From ( pq ) and 
tineutrinos) is 



= F^ + F^ + Fi, + F, 

= (1-p)F° + (2+p + p)F°+(1-p)F0 



Summarizing, the equations 



59| , |63D can be written in the compact notation 

(64) 



f p 1 — p 
p 1 — p 
\l—p 1—p 2+p + p 




The above general expression holds for both, the normal and the inverted mass hier- 
archy. The survival probabilities p and p for the normal hierarchy are given in (|55D 
and (|60|) respectively. These probabilities depend on the specific parameters of the 
masses and mixing scheme. 
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3.3 Earth matter effects on the u^. spectrum 

The neutrino trajectory inside the earth before reaching the detector depends on 
the direction of the supernova relative to the earth and the time of the day. The 
comparison of signals from different detectors would allows one to reveal the earth 
matter effects. Also certain features of the energy spectra can reveal the earth matter 
effect even from the observations in one detector. 

The mass eigenstates arriving at the surface of the earth oscillate in the earth 
matter. Let Pie be the probability that a mass eigenstate z/j entering the earth reaches 
the detector The flux of z/e at the detector is 

Fj^ = Y.P^eF^ ■ (65) 

i 

Inserting Fj from (0), we get 



where a^'s are as defined in (|52D , and we have used the unitarity condition Pie = 1- 



Thus, the Uf. flux at the detector can be written as 

FF=p''f!+{1-p^)F'^ , (67) 

with 

p'' = Y.a,P,e . (68) 

i 

Comparing (67) with F^ (^41), we find the difference in the Ve fluxes at the detector 
due to the propagation in earth equals 

FF-Fe = {p''-pm-F'^) , (69) 

where p is given in ([55|) . The earth matter effect can be quantified by the difference 
of probabililty {p^ — p) : 

P"" -p = Y.('^iP^e-\Ue^?) ■ (70) 



Using definitions in (|52D and YliPie = 1, we can write explicitly 

p" -p = PH{P2e-\Ue2?){l-2PL) + {Pze-\Ue?.\^){l-PH-PHPL) ■ (71) 



The second term in ( |7TD can be neglected. Indeed, inside the earth, z/3 oscillates with 
a very small depth: 



Pze-\UJ'< {'^-^^]s\r?29e, , (72) 

where Vearth is the effective potential of in the earth. For neutrino energies of 5 - 
50 MeV, we have 2EVearth/ ^rnl^^^ < IQ-"^. Moreover, sin2 2^e3 < 0.1, so that 

P3e-|t/e3|'<10~' • (73) 
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Finally, we can write 



P 



D 



2e 



|f/e2p)(l-2Pz.) 



(74) 



In general, when the signals from two detectors Dl and D2 are compared, we get the 
difference of fluxes 



D2 



(75) 



where P^^J and p!^ are the ^ ^2 oscillation probabilities for the detectors Dl and 
D2 correspondingly. 

According to (^) the earth matter effect is factorized: it is proportional to the 
difference in the original Ve and Vx fluxes, the conversion factor Ph{^ — ^Pl) inside 
the star, and the difference of earth oscillation probabilities at the two detectors. 
Let us consider these factors separately. 

1) — P°: Since the z/g spectrum is softer than the spectrum, and the lu- 
minosities of both the spectra are similar in magnitude the term {Fj^ — F^) is 
positive at low energies and becomes negative at higher energies where the flux 
overwhelms the z/g flux. Therefore, the earth effect has a different sign for low and 
high energies, and there exists a critical energy E^, such that F^{Ec) = F^{Ec), where 
this change of sign takes place. Since the cross section of the neutrino interactions 
increases with energy, the Earth effect is expected to be more significant at higher 
energies (if all the other factors are only weakly sensitive to the neutrino energy). 

2) Pff(l — 2Pl): This factor characterizes the neutrino conversions inside the star. 
Ph can be looked upon as a suppression factor due to the conversions at the higher 
resonance. Indeed, in the limit of Ues and Ph 1, ([7^) reduces to the expression 
for the earth effects in the case of two neutrino mixing : 



[p'' -p]2u = {P2e-\Ue2n{l-2PL) 



(76) 



which is equivalent to the one used in literature |^0[ in the context of day-night effect 
for solar neutrinos. Therefore (|7^ can be looked upon as the earth matter effect 
due to the two neutrino mixing ( |7^ ) suppressed by a factor of Ph- The mixing of 
the third neutrino thus plays a major role, making the expected earth effects in the 
case of supernova neutrinos smaller than those expected in the case of solar neutrinos 
for the same mixing scheme and in the same energy range. If the if-resonance is 
completely adiabatic, the earth effect vanishes: all the z/g produced are converted 
to z/3 in the star, and the earth matter effect on 1/3 is negligibly small (as we have 
established through (|72D ). 

3) -PiP ~ -fie''- If thfi neutrino trajectory crosses only the mantle of the earth, one 
can use a constant density approximation which gives 



p(i) _ p 

2e 2e 



(2) 



sin 2^^ sm{292-29e2) 



sm 



ndi 



— sm 



(77) 



Here 6*^ and I 



are the mixing angle and oscillation length inside the earth respec- 
tively, and di is the distance travelled by the neutrinos inside the earth before reaching 
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the detector Di. The first two terms on the right hand side of ( [77D are positive for 
the scenarios with the SMA and LMA solutions, so that the sign of (P2e'' — P^) is 
the same as the sign of the term inside the square bracket in (|77D . For the scenario 
with the VO solution, the earth matter effects are negligible: for Am? ~ 10~^° eV^, 
the mixing in the earth matter is highly suppressed, so that sin 2^^ in (|77| ) is very 
small. 

If neutrinos cross both the mantle and the core, the parametric enhancement of 
oscillations may occur, which leads to the appearance of parametric peaks apart from 
the peaks due to the MSW resonances in the core and the mantle [0]. Correspond- 
ingly the factor {P2^f} — P^) will be a more complicated function of the neutrino 
energy. 

To summarize, the earth matter effects on the z/g spectrum can be significant 
only for the scenarios with the SMA or LMA or LOW solutions. Moreover, the H 
resonance needs to be non-adiabatic in the case of normal mass hierarchy. 



3.4 Earth matter effects on the z/g spectrum 

The oscillation effects, are determined by the survival probability of the electron 
antineutrinos, , at the detector. Considerations similar to those for ([60| ) lead to 



Pleil-Pl 



P2ePL 



Then from Eqs. (|59|), (|60|) and (^), we obtain 



F^-F,= (Pie - M) ■ (1 - 2P,) . (F° - F°) 



(7J 



(79) 



where we have neglected the oscillations of z/3 inside the earth. Generalizing the result 



m 



we find the difference in the fluxes at two detectors Dl and D2: 



e e 



D2 



iHl^ - Hf ) ■ (1 - 2Pl) ■ (F° - F° 



^0) 



The earth matter effects for the antineutrinos are thus also factorized: they are 
proportional to the difference in the original z/g and fluxes, the factor of (1 — 2Pl) 
which depends on the conversions inside the star, and the difference [Pil^ — P^) 
between the oscillation probabilities inside the Earth for the neutrinos reaching the 
two detectors. Note that due to the absence of the i7-resonance in the antineutrino 
channel, there is no suppression factor similar to Ph ( |75D in the neutrino case. 

The factor [F^ — F^) is positive at low energies and negative at high energies. 
The matter effects then change sign at an energy Ec where Fg{Ec) = F^{E^. Since 
the neutrino cross section increases with energy, the observed effect is expected to 
be larger at higher energies (if all the other factors are only weakly sensitive to the 
antineutrino energy). 

In the approximation of a constant density, the factor depending on the oscillation 
probabilities inside the earth is 



P'S - Pie ~ - sin 20,^ sin(2^S - 26^2] 



,(2) 



sm 



sm 



^1) 
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where 6^ is the mixing angle inside the earth for the antineutrinos. For the an- 
tineutrino channel 9^ < ^ 1 for SMA solution and 9^ is strongly suppressed by 
matter in the VO case. Therefore the earth matter effects on the Ue spectrum can 
be significant only for the scenario with the LMA (as well as LOW) solution. In this 
scenario, sin 2^^ > and sin(2^^ — 29e2) < 0, so that according to (0), the sign of 
{Pil^ — -Pie'') is the same as the sign of the oscillation term inside the square bracket 



in Bl 



To summarize, in the case of normal mass hierarchy, the earth matter effects on 
the Ug spectrum are significant only for the scenario with the LMA solution. The 
effects are practically unaffected by the mixing of z/3 and change sign at an energy Ec 
such that FP(Ee) = F°(^c)- 



4 Effects of neutrino conversion for the mass spec- 
tra with normal hierarchy 

In this section, we shall consider the neutrino conversion effects for specific 3z/ schemes 
with the normal mass hierarchy. The general expressions for the transition proba- 
bilities and the neutrino fluxes at the detectors are given in Sec. |^, where they are 
expressed in terms of the total survival probabilities p and p of electron neutrinos 
and electron antineutrinos respectively. The values of p and p need to be calculated 
separately for each specific scheme. Notice that, given a scheme, the least known 
parameter is |t/e3p. We shall consider the effects of neutrino conversion in the three 
possible regions of |f/e3p: these correspond to the regions I, II and III as described in 
Sec. (see also Figs. ^ and ||.) 



4.1 The scheme with the SMA solution 

The mass and flavor spectrum of the scheme is shown in Fig. ^. The non-electron 
neutrinos and Ur mix strongly in the mass eigenstates z/2 and z/3. The electron 
flavor is weakly mixed: it is mainly in ui with small admixtures in the heavy states. 
The solar neutrino data is explained via the small mixing angle MSW solution, with 
the parameters as given in (^). The level crossing scheme is shown in Fig. |^a. Both 
the resonances are in the neutrino channel. 
In this scheme, 

|t/ei|'~l , |t/e2r~^sin2 2^0~lO-3 . (82) 

Let us first consider the antineutrino channels. There is no resonance here, and 
the mixing in matter is suppressed. As a consequence, the adiabaticity condition is 
fulfilled in the L layer and ~ 0. Then according to the level crossing scheme 
(Fig. ^), the following transitions occur inside the star: 

I^e ^ Z/i , Z/2 , Ur' ^ Z>3 • 
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Figure 6: Neutrino mass and mixing pattern for the scheme with the small mixing 
angle MSW solution of the solar neutrino problem. The boxes correspond to the mass 
eigenstates. The sizes of different regions in the boxes show admixtures of different 
flavors. Weakly hatched regions correspond to the electron flavor, strongly hatched 
regions depict the muon flavor, black regions represent the tau flavor. 



Using |?7e2p 1 and ~ 0, we get from (|60D the survival probability for z/, 



P 



\UelWl-PL)-l 



^3) 



Let us consider now neutrino channels. The value of the survival probability p 
depends on the region in which the oscillation parameters of the H resonance lie. 

Region I: Since the H resonance is adiabatic, the level crossing scheme (Fig. ^) 
leads to the following transitions: 



T^e 1^3 

At Ph ^0, the Eq. (H) gives 



P ~ If/esT < 0.03 



The flavor transitions are then complete, and independent of the adiabaticity at the 
L resonance. 

The final neutrino (and antineutrino) fluxes are 



The neutrino spectra have the following features: 

• The neutronization z/g peak disappears (it is suppressed by a factor of |f/, 



e3| 



^5) 



< 



0.03). Instead one would expect a neutronization peak which may be detected 
by the neutral current interactions. 
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• The antineutrino signal is practically unchanged. 

• The has the hard spectrum of F^, and therefore 

which is a clear signal of mixing. 

• The z/j. spectrum is composite, i.e. it contains both the soft (original z/g) and 
the hard (original Vx) components. 

The earth matter effects on the z/g spectrum (|75D are suppressed by the factor of 
Ph ~ 0. Also, since the mixing in the antineutrino channel is suppressed in matter, 
the earth matter effects on the Ve spectrum are negligible. As a result, one expects 
practically the same neutrino (and antineutrino) signal in all detectors. 

Region II. Neutrinos jump partially between the third and the second levels in 
H region, so that the following transitions occur: 

Z/e Z/i, Z/2, Z/3 , 1^1, 1^2 , l^r' 1^1, 1^2, 1^3 ■ 

The z/g survival probability p is given by the general formula (|55D. The value of Ph, 
and hence the survival probability p, depends on the neutrino energy, although this 
dependence is relatively weak, as demonstrated in Figs. and 0. 
The features of the spectra are: 

• the neutronization peak is distributed in both, z/g and z/^,, 

• the z/g and z/^, spectra are composite, 

• the z/g spectrum is unchanged. 

The earth matter effects can be observed in the z/g spectrum. The magnitude of the 
effects is proportional to Ph ([75|). For z/g spectrum the effect is negligible. 

Region III. No conversion occurs at the H resonance, so that in the star 

z/g Z/i, Z/2 , Z/^/ Z/i, Z/2 , Ur' ^ ■ 

From (^) and (|82|) we get using Ph ~ 1: 

P^\Uel\'PL+\Ue2\\l-PL)^PL • (86) 

The features of the spectra are the same as those described above in the case of region 
II. 

The earth matter effects on the z/g spectrum in the SMA scheme can be significant 
in regions II and III, where the factor of Ph is of the order of unity. Indeed, in ([77[), 
the value of (P2e — |f^e2p) can be as large as 0.25 in the energy range of 20-40 MeV, 
where the term {F^ — F^) is also significant. 

In Fig. 1^, we show the z/g spectrum for different distances travelled by the neutrinos 
through the earth. As follows from Fig. ^ the net effect is < 10% even with optimistic 
values of Ph, Pl and 6'g2. 
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Figure 7: The earth matter effects on the z/e spectrum for Ph = 1 and Pl = Qin the 
scheme with the SMA solution (Am^ = 10"^ eV^ , sin^ 26'0 = 10"^). The dotted, 
dashed and sohd hues show the number oi v^, — N charged current events when the 
distance travelled by the neutrinos through the mantle of the earth is ci = km, 
d = 6000 km, and d = 10000 km respectively. 

4.2 The scheme with the LMA solution 

The mass and flavor spectrum of the scheme is shown in Fig. |^. The solar neutrino 
data is explained via the 1^2 resonant conversion inside the sun with a large 

vacuum mixing angle (^). The level crossing scheme is shown in Fig. It differs from 
the previous one by the mixing at the L-resonance. Since the vacuum mixing angle 
^0 is large, there is a significant Ue ^ 1^2 conversion even though the antineutrinos do 
not encounter any resonances. The evolution is adiabatic in the antineutrino channel 
in both the H and L layers. The following transitions occur: 



1^1 



1^2 



Using ~ 0, the z/g survival probability (|60|) equals 



el 



COS Oct 







^^3 



^7) 



Let us consider transitions in the neutrino channels. The L-resonance is adiabatic 
(Fig. Hp, Pl ~ 0, and from (pH), we get the z/g- survival probability 



P 



\Ue2VPH+\Ue,\\l-PH) 



sin^ OqPh + W, 



e3| 



H 



Notice that, depending on the value of Ph-, the Vf.- survival probability takes the 
values between 

(89) 



\UeZ? <P<\Ue2? 
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Figure 8: Neutrino mass and mixing pattern of the scheme with the large mixing 
angle MSW solution for the solar neutrino problem. 



The fluxes at the earth are determined by (^), with p and p given in (^) and 
(87). Their features depend on the adiabaticity at the H resonance, which is decided 
by the region in which |f/e3p lies. 

Region I. The level crossing scheme (Fig. ^) with adiabatic transitions in the 
H resonance layer leads to the following transitions inside the star: 

Z/g ^ 1/3 , Z/i , U-r' Z/2 . 

Since Ph ~ we get p ~ jt^esP, and the final spectra at the earth: 

Fg ^ cos' ^oF° + sin^ ^0F° , 

^ F° + sin' ^qF" + (2 + cos' ^0)^0 . (90) 

The following features can be observed: 

• The neutronization peak almost disappears from the channel (suppressed by 
a factor of jt/esl' ^ 0.03) and appears in the Ux channel. 

• The Ue spectrum is composite: A mixture of the original spectrum and the 
original harder spectrum, is determined by the solar neutrino mixing angle. 
This feature distinguishes the LMA scheme from the one with the SMA solution. 

• The z/g spectrum is hard, practically coinciding with the original spectrum. 
This leads to 
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• The spectrum contains components of all the three original spectra {ve-, ^e-, ^x)- 

No earth matter effects are expected in the z/g spectrum, since Pu ~ 0. At the same 
time, the z/g spectrum can show significant earth matter effects. This is an important 
signature of the scenario. 

Region II. The jump probability in H resonance is substantial and the following 
transitions occur: 



The final spectra are determined by (|^) with p given by (^) : 

p ^ sin^ 6q Ph ■ 

Depending on Pjj and sin^6'0 ~ |f^e2p ~ 0.2 — 0.4, the probability p can take values 
between 0.02 and 0.4. 

The final spectra can be written as 



^ sm'e^PHF^ + il-sm^eMF: 



F-e ^ cos^ ^0F° + sin^ ^0F° 
4F, ^ (l-sin2^0Pj,)F° + sin2^0F° + (3-sin2^0 + sin2^0Pj,)F° . (91) 

The features of the spectra are: 

• The neutronization peak consists of both, v^, and fluxes, with the proportions 
determined by 9q and Ph ■ 

• The z/e and z/g spectra are composite, with admixture of the hard component 
determined by the solar mixing angle. 

• The spectrum contains components of all the three original spectra (z^e, i^e, ^x)- 

The earth matter effects can be significant for both the z/g and z/g spectra. The effects 
on the z^e spectrum are suppressed moderately by a factor of Ph- 



Region III: Here Ph ~ 1, so that according to (|88D , p ^ |?7'e2p ~ 0.2 — 0.4. The 



if-resonance is inoperative, and the following transitions occur: 

z/e ^ Z/2 , Z/^' Z/1 , Ur' ^ Z/3 . 

The features of the spectra are the same as those described in the case of the region 
II. Since Ph ~ 1, the earth matter effects are expected to be significant for both z/g 
and spectra. 



Let us consider in some details the the earth matter effects which can give impor- 
tant signature of the scheme under consideration. For the z/g spectrum the effect can 
be significant in regions II and III. Inserting = in (|77D , we get 

- Pe ^ Ph ■ (P2e - \Ue2\') " (F° - F°) (92) 
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Using the constant density earth approximation ( |T^ ) we estimate that the factor 
(P2e ~ I^e2p) Can be as large as 0.3 in the energy range of 20 - 50 MeV. 
Similarly, for the antineutrinos, we get 

Fi^-F,^(Ae-|f/eir)-(F°-F.°) . (93) 

In the antineutrino channel, the effect can be significant in all three regions of f/es, 
since there is no Ph suppression, and again the factor (Pie — If/eiP) can be as large 
as 0.3 for energies around 20-50 MeV. 



N N 




(a) (b) 

Figure 9: The earth matter effects on (a) the i^e spectrum and (b) the spectrum for 
= 1 in the scheme with the LMA solution (Am^ = 2 ■ 10"^ eV^ sin^ 26^0 = 0.9). 
The dotted, dashed and sohd lines show the spectra of the number of u — N charged 
current events when the distance travelled by the neutrinos through the earth is ci = 
km, d = 4000 km, and d = 6000 km respectively. 



In Fig. ^ we show the u^, and z/g spectra for different distances travelled by the 
neutrinos through the earth. The sign of the earth effect (^, ^3|) is determined by the 
sign of (Pg° — P°) in the neutrino channel and (P° — P°) in the antineutrino channel. 
The effect changes sign at E = {E = E^ for the neutrinos (antineutrinos), as 
described in Sec. p| (Sec. 0). Since Pge > \Ue2? and Pie > |t/eiP, (|9|) and (H) 
imply that the earth effect results in the enhancement of neutrino (antineutrino) signal 
for E < Ec {E < Ec) and the depletion of the signal for E > Ec {E > Ec). Since 
the cross section of the neutrino detection interactions is larger at higher energies, a 
significant reduction in the Ue (i^e) flux is expected at energies greater than Ec (Ec). 
It may manifest itself as a dip in the final spectrum (see Fig. |]), and this spectral 
distortion can indicate the presence of significant earth matter effects, even with the 
spectrum observed in only one detector. 
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4.3 The scheme with the VO solution 



The mass and flavor spectrum for tfie scfieme is sliown in Fig. |T0[ Tlie solar neutrino 
data is explained via tlie vacuum oscillations with Am^ ~ 10~^° eV^ and a large 
mixing between and a combination of and (|^). 
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Figure 10: Neutrino mass and mixing pattern of the scheme with the vacuum 
osacillation solution for the solar neutrino data. 

The level crossing scheme is shown in Fig. |^c, where the resonance L is close to 
p = 0. It may occur in the neutrino or antineutrino channel. Depending on the details 
of the density profile in the star at low densities, it lies either in the non-adiabatic 
region (so that P^, ~ 1) or in the transition region. 

Let us first consider the antineutrino channel. According to the level crossing 
scheme, we get the transitions 

Ve Z/i, Z/2 , Vy^i Z/l, Z/2 , Vr' Vz ■ 

The final fluxes are given by (|55|) and for p we get from ( |60|) 

p = cos2^0(l-Pi) + sin2^0Pi . (94) 

This implies that the value of p lies between cos^ and sin^ 6'©, and the z/g spectrum 
is composite. 

The fluxes in the neutrino channel depend on the region in which the oscillation 
parameters of the H resonance lie. 
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Region I. The H resonance is adiabatic, and the following transitions occur in 
the neutrino channel: 

Z/e ^ Z/3 , I/^/ Ui, U2 , Vr' ^ ^^1, 1^2 ■ 

Since Ph ~ 0, from (^) we get 

P ~ \Uez\^ < 0.03 . 

The Vf, then has the hard spectrum of the original v^. The observed features of the 
final fluxes are the same as those in the scheme with the LMA solution and an adi- 
abatic H resonance (|90|). The earth matter effects are small for both the Ve and z/g 
spectra, since the mixing angle 9e2 is suppressed in the earth matter. 

Region II. Jump probability in the H resonance is substantial, so that the tran- 
sitions inside the star are: 

Ve 1^1, 1^2, 1^3 , l^fi' 1^1, 1^2 , l^T' 1^1, 1^2, 1^3 ■ 

From (^), we can infer 

|f/e3r <P< COS'^0 . (95) 

where the lower bound corresponds to Ph = 0, whereas the upper bound is for 
Ph = Pl = 1. The z/g spectrum is then composite. 
The features of the spectra are: 

• The neutronization peak should be observed both in z/g and in channels. 

• The Ue spectrum is composite, the admixture of the original hard component is 
determined by the solar mixing angle as well as by the degree of adiabaticity in 
the H and L conversion regions. 

• The Ue spectrum is composite, the admixture of the original hard component 
is determined by the solar mixing angle as well as the adiabaticity in the L 
resonance region. 

• The I'x spectrum contains components of all the three original spectra (z/g, i^e, ^x)- 
The earth matter effects are small for both the z/g and the z/g spectra. 

Region III. Since the H resonance is completely non-adiabatic, the neutrino 
transitions take place in L resonance region only: 

Z/g Z/i, Z/2 , Z/^' Z/i, Z/2 , Vr' ^ ^^3 • 

From (^) we get for Ph ~ 1: 

P ~ \Uei?PL + \Ue2Wl - Pl) ~ cos^ BqPl + sin^ ^0(1 - Pl) . (96) 
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Thus, 

sin^ 9q<p< cos^ 9q . (97) 

The i>e - spectrum is composite. The features of the spectra are the same as those 
observed above in region II. The earth matter effects are small for both the z/g and 
the z/g spectra. 



The bi-maximal mixing scheme 




Figure 11: The neutrino mass and mixing pattern of the "strict" bimaximal mixing 
scheme. 

An extreme case of the vacuum oscillations is the bimaximal mixing scenario 



The mass and flavor spectrum is as shown in Fig. [T^. and z/^ mix maximally in 

{v^ — i't)/ \/2 maximally 



/ -\/2; in turn, the orthogonal combination z/^ 



mixes maximally with z/g. In this scenario, Uez = 0. The mass eigenstates are then 



J^i = [J^e - v)/v2 , Z/2 
The following transitions take place: 

Z/e^Z/i,Z/2 , Z/^' 

Z/g " " 



1^3 



Z/i, Z/2 , Ur' —> 1^3 , 
Z/i,Z/2 , Ur' 1^3 • 

The mass eigenstates z/i and z/2 (z/i and z/2) arriving at the earth are observed with 
equal probability as electron or non-electron neutrinos. The expected fluxes are then 

Fg = (l/2)Fg° + (1/2)F: , 
F, = (1/2)F°+(1/2)F° 



4F. 



(l/2)Fg° + 3F,° + (l/2)F° 



(98) 
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The neutrino spectra are affected in the following way: 

• The neutronization peak consists of equal fractions of electron and non-electron 
neutrinos. 

• The z/g-spectrum is "composite" , it has both the soft and the hard components 
in equal proportions. 

• The z/g-spectrum is also composite with both the soft and the hard compo- 
nents in equal portions. This spectrum has a lower average energy than the z/g 
spectrum, but the high energy tails of both the spectra are identical. 

• The spectrum has a mixture of all the three initial spectra. 

The earth matter effects are small for both z/g and z/g spectra due to the small value 
of Am| ~ 10-1° eV2. 

5 Conversion probabilities and neutrino fluxes with 
the inverted mass hierarchy 

In the case of the inverted mass hierarchy ([T9|), the mass eigenstates vi and V2 are 
heavy and degenerate, whereas the third state z/3 is much lighter: mi ~ m2 ^ m^. 
(see e.g. Fig. [T^. for the scenario with the SMA solution). 

As shown in Figs. §b and |^d, the if-resonance is in the antineutrino channel, 
whereas the L-resonance lies in the neutrino channel for the SMA and LMA solutions. 
In the case of the VO solution, the L resonance can be in either the neutrino or the 
antineutrino channel. 

Like in the case of the normal mass hierarchy, the fluxes at the earth can be 
written in terms of the survival probabilities p and p of z/g and z/g respectively, as 
given in ( |64D . Let us find the expressions for the survival probabilities p and p. 

As follows from the level crossing scheme (Figs. §a and|^d), the matter eigenstates 
in the high density region (p ^ Ph,Pl) are 

^Im = Z^m' ' ^2m = , J^3m = ^r' , (99) 
l^lm = i^T' , l^2m = t'/x' , i^3m = ■ (100) 

The state z/g, which coincides with in the production region, crosses the H- 
resonance layer adiabatically (since the resonance is in the antineutrino channel and 
the mixing in the neutrino channel is small), and reaches the L-layer as v^m- It reaches 
the surface of the earth as z/i with the probability and as with the probability 
(1 - Pi). Then 

P=|t/el|'PL+|f/e2|'(l-PL) • (101) 

The state z/g, which coincides with z/sm in the production region, flips to z/i^ in 
the resonance region H with the probability P^. Depending on the parameters of the 
solar neutrino solution, the propagation of the antineutrinos in the L-resonance region 
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Figure 12: The neutrino mass and mixing pattern for the scenario with inverted 
mass hierarchy and the SMA solution for the solar neutrino problem. 



{i. e. near p = 0) may or may not be adiabatic which is described by P^. Computing 
the probabilities of transitions of to ui, z/2, and performing projection of the 
result back onto back we get (similarly to (|55[) the survival probability p for the 
z/e: 

P = \Uel\^PH{l - Pl) + \Ue2?PHPL + WezW^ ~ Pr) • (102) 

(Here we have also averaged out the interference effects between the mass eigenstates, 
so that the three terms in (|1U2D correspond to the three mass eigenstates.) 



Let us consider the earth matter effects on the z/g and Vf. spectra for the inverted 
mass hierarchy. Since there is no if- resonance in the neutrino channel, the earth 
effects for the inverted hierarchy are the same as those described for the normal mass 
hierarchy (see Sec. |3.3|) with Ph = 1: 



The matter effects are significant for the scheme with the LMA solution. The sup- 
pression factor of Ph [that was present for the normal mass hierarchy - see eq. ([73D] 
is now absent, so the earth matter effects can be larger than in the case of the normal 
mass hierarchy (see Fig. ^a). 

The earth matter effects on the antineutrino spectra can be calculated as follows. 
Using the same arguments as for (|102|) , we find the survival probability of Vf. at the 
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detector: 

= PlePnil - Pl) + P2ePHPL + P?M " Pr) • (104) 

Then (|0|) and (|10|) give 

- P = (Pie - |f/ein(l - 2Pl)P// + {he ' \Ue^?){l ' Pr - PrPl) , (105) 

where we have used Pie = 1- Since z/3 oscillates inside the earth with a very small 
depth [inequality (|7^) is valid with P^e replaced by P^e] , the second term in ( |105| ) can 
be neglected. Therefore, finally we get 

pm _ pm ^ p^^pii) p(2)) . (1 _ 2p^) . ^po _ po^ ^^gg) 

The earth matter effects on the Ue spectrum for the inverted hierarchy are then the 
same as those described for the normal hierarchy in Sec. |3.3| , but further suppressed 
by a factor of Ph- Significant effects can be observed only with the LMA scenario 
and when the H resonance is non-adiabatic. 

The earth matter effects on both and z/g spectra are factorized, as in the case of 
the normal mass hierarchy. They are proportional to (i) the difference in original z/g 
(z/g) and fluxes, (ii) a factor which determines the conversion inside the star, and 
(iii) the difference of the oscillation probabilities inside the earth. 

Notice that the expressions for the earth matter effects in inverted mass hierarchy 
106|) have the same form as the expressions in the normal hierarchy case 



case 



^ with substitution: z/ z/, Ph ^ Ph- 
Let us now calculate the values of the survival probabilities p and p for the three 
scenarios with the SMA, LMA and VO solutions for the solar neutrino problem re- 
spectively. Within each scenario, the flnal spectra would depend on the adiabaticity 
of antineutrino conversions in the H resonance layer, i.e. on the value of |f/e3p- 

We shall consider the effects of neutrino conversions in three possible regions of 
|f/e3p: regions I, II and III as described in Sec. ^.5|. When |t/e3p lies in region III i.e. 



the H resonance is inoperative, the resulting flnal spectra for a given solar neutrino 



solution are the same as those obtained in Sec. [4.1| , [4.2| and [4.3| for the normal mass 
hierarchy with the same mixing parameters. Therefore, we need to consider only the 
cases when |t/e3p lies in regions I or II. 

5.1 The scheme with the SMA solution 

The characteristics of the z/g spectrum are independent of the conversions at the H 
resonance layer [see eq. (|1U1|)]. The L resonance is in the transition region so that 



partial jumps from one level to another occur. As a consequence, from the level 
crossing scheme (Fig. ^d), we get the following transitions: 



Z/e Z/i, Z/2 , Z/^' Z/i, Z/2 , U^' ^ Z/3 . 

|2 



Neglecting terms proportional to the small |f/e2| ( |101[ ), we get 

Pl ■ 
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The Ve spectrum is thus composite. 

Let us now examine the final v^. spectrum. Since the L resonance is in the neutrino 
channel, the antineutrino transitions in the L resonance layer are adiabatic, so that 
Pl ~ 0. The characteristics of the final z/g spectrum depend on the region in which 
the value of |C/e3p lies. 

Region I. The following transitions take place inside the star: 

Ve-^yZ , V^, Z/2 , Vr' ^ V\ ■ 

For Pi^ ^ the Eq. ( |T02D gives 

V ~ |f/e3p < 0.03 . 

The final v^, spectrum is then almost completely the original hard spectrum of v^: 

F- ^ F° 

The features of the final spectra are: 

• The neutronization peak contains both, Vf. and fiuxes. 

• The Vf. and v^. spectrua are composite. 

• The Vf. spectrum is practically the original hard spectrum of v^. 



The earth matter effects on the z/g spectrum ( p.03| ) are small: < 10% (see Sec. [4.1[ ) 



but still may be observable. The earth matter effects on the spectrum ( |106| ) are 



negligible. They are suppressed by the factor of Ph ~ and further by the small 
mixing of the antineutrinos in the earth matter. 

Region II. Taking into account jumps of states at H resonance we get the fol- 
lowing transitions: 

Ve Z/i, Z/3 , V^, z/2 , Ur' 1^1, 1^3 ■ 

Inserting in eq. ( |102| ), we find 

Ph cos^ 6q ^ Ph . 
The following features are then expected in the final spectra: 

• The neutronization peak contains both Uf. and fiuxes. 

• All the three spectra, of z^e, I'e and z/^., are composite. In contrast with previous 
case (region I) now z/g spectrum is also composite. 



The earth matter effects (|103| ) are small: < 10% (see Sec. [4.1| ) for z/g, and negligible 
for the Up. 
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5.2 The scheme with the LMA solution 



The neutrino transitions at the L resonance are adiabatic, so that the following tran- 
sitions take place according to the level crossing scheme (Fig. ^jd): 

Inserting ^ in Eq. ([101|), we get 

P ^ |t/e2r ~ sin^ ^0 ~ 0.2 - 0.4 , 

and consequently, the spectrum is composite: 

F, ^ sin^ OqF^ + cos^ OqF^ . 

Let us examine the final spectrum. The antineutrino transitions in the L res- 
onance region are also adiabatic, so that Pl ~ 0. Conversion in the H resonance 
region, and therefore the characteristics of the final spectrum depend on |t/e3p- 

Region I. Here the H resonance is adiabatic and we get the following transitions: 

Since Ph ^ 0, the Eq. ( |T0^ ) gives 

P ~ It^esT < 0.03 . 

and consequently, the Ue spectrum is then almost completely the original hard spec- 
trum of z/^,: 

F- ^ F° 

The features of the final spectra are the same as those of the scheme with the inverted 
mass hierarchy and SMA solution (with |?7e3p in region I). The earth matter effects 
on the z/e spectrum (|103|) can be significant, as can be seen from Fig. |^a. In contrast, 
the earth effects on the z/g spectrum ( |106| ) are negligible: they are suppressed by the 
factor of Ph ~ and by the small mixing of the antineutrinos in the earth matter. 

Region II. The transitions in the H resonance layer are incomplete and therefore 
according to the level crossing scheme we have 

z/e Z/i, Z/3 , Z/^/ Z/2 , Z/^' ^ Ui, . 

Inserting ^ in eq. (|102|) , we get 

PnlUeil^ = Pncos^ee . 

The features of the final spectra are the same as those of the scheme with the inverted 
mass hierarchy and SMA solution (with |t/e3p in region II). The earth matter effects 
can be significant both in the z/e ( [103| ) and in the z/e channels ( |106| ) (see Fig. |^). 
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5.3 The scheme with the VO solution 

The conversions in the H resonance layer does not influence z/g flux. In the L resonance 
layer the adiabaticity is moderately broken: Pl ^ 0. From the level crossing scheme 
(Fig. the transitions taking place inside the star are: 

Z/e Z/i, Z/2 , Z/^' Z/i, Z/2 , U^' ^ • 

From ( |101| ), we get 

sin^ 9q < p < cos^ 6*0 . 

and since sin^6'0 > 0.3 we get p = 0.3 — 0.7. The Ue spectrum is thus composite. 
Let us now examine the final Ue spectrum. 

Region I. Due to adiabaticity in the H resonance layer, the following transitions 
take place: 

Z>e ^ Z>3 , i^f,' i^l, ^2 , ^ Z>1, V2 ■ 

Here Ph ~ 0, and from ( p.02|) we get 

It/esr <0.03 ■ 

The z^e spectrum is then almost completely the original hard spectrum of z/^.: 

F- 

The features of the final spectra are similar to those of the scheme with the inverted 
mass hierarchy and SMA solution (with |f/e3p in region I): the spectrum is com- 
posite whereas the z/g spectrum is purely hard. The earth matter effects on both z/g 
and spectra are negligible. 

Region II. The adiabaticity is moderately broken both in the H and in the L 

resonance layers and the transitions proceed as 

Ve Z^2, , i^l, 1^2 , l^r' ~^ 1^1, 1^2, 1^3 ■ 

Correspondingly the Eq. (|102|) gives 

p ^ Ph[cos^ 9^(1 - Pl) + sin^ O^Pl] . 

The features of the final spectra are the same as those of the scheme with the inverted 
mass hierarchy and SMA solution (with |f/e3p in region II). The earth matter effects 
on both z/e and Ve spectra are negligible. 

Let us summarize the results. Some salient features of the final spectra with the 
inverted mass hierarchy are the following: 

• The neutronization peak contains both v^. and fluxes, in all the scenarios. 

• Also the final z/g spectrum is composite in all scenarios. The characteristics of 
this spectrum are independent of the oscillation parameters in the H resonance 
region. 
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• The final Vx spectrum is composite in all scenarios. 

• When the value of \Uez\^ lies in region I, the final Ve spectrum is almost the 
original hard spectrum of v^- For |t/e3p in regions II and III the z/g spectrum is 
composite. 

• The earth matter effects on the z/g spectrum can be significant for the scheme 
with the LMA solution, marginally observable for the scheme with the SMA 
solution and negligible for the VO solution. 

• The earth matter effects on the z/g spectrum can be significant only for the 
scheme with the LMA solution, and when the antineutrino transitions at the H 
resonance are not adiabatic. In all the other scenarios, the earth matter effects 
on z/e spectrum are negligible. 

• When the value of |t/e3p lies in region III, the final spectra (including the earth 
matter effects on them) are indistinguishable from those obtained with the 
scheme having the same oscillation parameters but the normal mass hierarchy. 



6 Signals of mixing and signatures of mixing schemes 



The future detection of a supernova neutrino burst by the underground neutrino 
detectors has been discussed in |5^, For a typical supernova at 10 kpc. 



about 5000 Ve events are expected to be detected at SK, and a few hundred events each 

and MACRO 



in SNO ) LVD ||58| and MACRO ||59|. These detectors can reconstruct the energy 
of the outgoing charged lepton from the Ve— P charged current interactions. SK and 
SNO can also observe the direction of the charged lepton. In addition, SNO can detect 
z/e and reconstruct its energy through the I'e — d charged current interaction inside the 
D2O. The feasibility of measuring the absolute values of the neutrino masses through 



the time-of-fiight delays has been studied in [O, IBO, p^. Here, we concentrate on 



the features of the final neutrino spectra that are relevant for identification of the 
neutrino mass and flavor spectrum. 

The effects of neutrino transitions on the final neutrino spectra can be observed 
through, e.g. (i) the partial or complete disappearance of the z/g neutronization peak, 
(ii) the apperance of hard / composite spectra of z/g or z/g, (iii) earth matter effects. 
In this section, we shall elaborate on the effects that can be observed in the earth 
detectors. 



6.1 The neutronization peak 

The neutronization peak can be identified as the burst in the neutrino signal during 
the first few milliseconds. If there is no neutrino mixing, the neutrinos will be pre- 
dominantly z/g. Let the expected ratio of the number of charged current (CC) events 
to the number of neutral current (NC) events during the burst be 

^0 - j^NC ■ 
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If neutrino conversions transform some Ve into v^i the number of charged current 
events decreases while the number of neutral current events remains unchanged. 
Then, the observed value of this ratio is < Rq . 

When the value of |f/e3p lies in the region I in the case of the normal mass 
hierarchy (so that the neutrino transitions in the H resonance layer are adiabatic), 
the z/e flux, and hence the number of charged current events, gets suppressed by a 
factor of ^ |t/e3p (see Sees. |4^ , 



^ ~ |?7e3|' < 0.03 . (107) 

The observation of Rb > 0.03 Rq would then point against |t/e3p in region I for the 
normal mass hierarchy. According to (^0]) , this will allow us to put an upper bound 
on the mixing parameter |?7e3p: 

It/esP < 3 ■ 10"^ . (108) 



6.2 The {E) inequalities and high energy "tails" 

In the absence of mixing, one should observe {E^J < (Ep^). The inversion of this 
inequality, i.e. {E^J > {EpJ is the signature of the i^-conversion. It implies that 
the contribution of the converted original hard u^. spectrum to the final z/g flux is 
signiflcantly larger than its contribution to the flnal Ue flux. 

In the case of the normal mass hierarchy, the inverted inequality for (E) will be 
observed when the transitions in the H resonance layer are adiabatic. The observation 
of the normal inequality {Ej^J < {EpJ can exclude this scenario. 

In the case of inverted mass hierarchy, the adiabatic transitions in the H reso- 
nance layer preserve the original inequalities for (E): {E^J < {EpJ. Therefore, an 
observation of the inequality {E,yJ > {EpJ rules out the scenarios which have the 
inverted mass hierarchy with the value of |f/e3p in region I. If the mass hierarchy is 
known to be the inverted one, an upper bound ( p,08|) on the mixing parameter |t/e3p 
can be obtained. 

The relative difference of the z/g, and spectra is especially signiflcant in the 
high energy ends ( "tails" ) of the spectra, where the fluxes decrease exponentially with 
the increase of energy. In the absence of mixing, one expects 

N,{E>Etaa) N,{E>Et,u) 

N,{E > Etail) ' N^{E > Etau) ^ ' 

for high enough energy Etaii (energy greater than the peak energies of the spectra). 
Here Ni{E > Etau) is the integrated number oi Vi — N charged current events above 
the energy Etau- 

The observation of the high energy ends of the z/g and Ve spectra, where the 
contributions due to the original z/g and fluxes would be negligible, allows one 
to measure the contribution of the original hard spectrum to the flnal z^g and z/g 
spectra. 
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Important conclusions can be drawn from the studies of the ratio 



R{Etail) 



(110) 




oo: 



This indicates that the tail of the Ue spectrum extends to significantly higher 
energies than that of the spectrum. This will testify for z/e ^ conversions. 
Moreover, it will be a clear signal that the Ve ^ conversions are negligible, i.e. 
the Ve spectrum is the original soft spectrum. This would be an indication of (i) the 
SMA solution in the conventional hierarchy (sec. [4.1|) , or (ii) highly non-adiabatic H 
transitions with the SMA solution in the inverted hierarchy (sec. p.l[ ). If the mass 
hierarchy is known to be the inverted one, then this signal indicates that the value of 
|[/e3p lies in region III. From (^l]), we obtain a very strong bound on the parameter 



Notice that the limit RiEtaii) —>■ would be true in the absense of any neutrino 
conversion at all, or in the absence of conversion of the electron neutrinos. This 
situation is not realized in any of the scenarios we have discussed. The observation 
of R{Etaii) would exclude all the scenarios under consideration and thus testify 
for the solution of the solar neutrino problem which differs from the MSW or VO 
solutions. Another possibility could be related to a compact star with a very sharp 
density profile, so that the L resonance is in the non-adiabatic region. In this case, 
the solar solution has to be SMA and the H resonance has to be inoperative for 
neutrinos. 

3) R{Etaii) const: 

This would testify for the contribution of the original z/^. to both final z/g and z/g 
spectra. Moreover, the value of the constant will give us information about the mix- 
ing parameters. 

In principle, the spectra of can be reconstructed by having a series of neutral 
current detectors with different thresholds. The spectrum will always have a 
dominating hard component. Therefore, the comparison of the high energy tails of 
the z/e, z/e with that of can directly give the measure of the hard component of the 
z/e, z/g spectra. The ability to measure the z/^. spectrum would also be useful in order 
to check for the presence of sterile neutrinos [0. 

6.3 The identification of a composite spectrum 

The final z/e (z/e) spectrum can be qualitatively divided into three types: (a) the 
original "soft" spectrum of the z/g (z/g) [corresponding to the survival probability p = 1 
{p = 1) which would be the case in the absence of conversion], (b) the "hard" original 
spectrum of [corresponding to the survival probability p = (p = 0) which would 




e3 




(111) 



2) R{Eta^l) 0: 
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be the case when there is a complete interchange of spectra], and (c) the "composite" 
spectrum, which is a mixture of the original soft and the hard spectra in comparable 
proportions. 

In order to distinguish between (a) and (b), the values of "temperatures" (or 
average energies) of the original spectra need to be known, or the and Ve spectra 



have to be compared (as in Sec. |6.2|) . A composite spectrum may be identified without 
the knowledge of the initial temperatures and independently of the observation of the 
other spectrum, through the "broadening" phenomenon. 



As described in Sec. the instantaneous original neutrino spectra have a nar- 
rower energy distribution than the Fermi-Dirac one (see Fig. ||). The mixing between 
two neutrino species with different mean energies gives rise to an effective broadening 
of the spectrum (see Fig. |T3|), so that The final spectrum need not be pinched even 
though the two individual original spectra were pinched. In other words, the effec- 
tive 7] for a mixed spectrum may be negative even though the individual ?7j's of the 
constituent spectra were positive. The broadening of the spectrum can be checked 
by fitting the final spectrum with the parameters Tj and r/j (|^) and establishing the 
sign of r]i. 




5 1 1.5 2 2.5 



E 



3 <E> 



Figure 13: The number oi v — N charged current events for (A) a Fermi-Dirac 
spectrum (T = 3,// = 0): the sohd line, (B) a pinched spectrum (T = 3,?] = 3): the 
dotted line, and (C) a mixture of two pinched spectra (T = 3, = 3 and T = 8, r/ = 1 
with Ena = 5 MeV): the dashed line. 

Note that the "pinching" phenomenon is established in a model independent man- 
ner only for a constant mean energy of the spectrum. A considerable variation in the 
average energy of neutrinos of a given species during the cooling phase could be re- 
sponsible for the broadening of the time integrated spectrum even in the absence of 
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mixing. Therefore, in order to establish the broadening due to mixing, it is crucial to 
study the spectra in short intervals of time. 

Clearly, it will be difficult to observe broadening of spectra, if the mean energies 
of the original spectra are very close, or if one of the original spectrum dominates in 
the final mixed spectrum. 

The observation of broadening of the z/g spectrum would point against an adiabatic 



H resonance in the normal hierarchy (see sees. |4.1| , [4.2| , |4.3| ). This would imply a 
strong upper bound (|108|) on the value of \Ue. 



^e3| 

The observation of broadening of the Ue spectrum would be a strong indication 
against 

• the SMA solution with the normal hierarchy (Sec. 



completely adiabatic H transitions in the case of inverted hierarchy (sees. ^T 
31), and 



highly non-adiabatic iJ-transitions for the SMA solution with the inverted hi- 
erarchy (Sec. [SH] ). 



In particular, if the mass hierarchy is inverted, the broadening of z/g spectrum indicates 
that the H resonance is not completely adiabatic. This gives a strong upper bound 
(|108D on |f/,3p. 



If both the Uf, and u^, spectra are established to be composite, then the upper 
bound (|108|) on |^7e3p holds irrespective of the type of mass hierarchy. 



6.4 Earth matter effects 

The earth matter effects can be observed through the comparison between the spectra 
at two detectors, or through the studies of the distortion of a spectrum in one detector 
(See Figs. 0, |). 

The earth matter effects are significant in the neutrino channel only in the cases 
of (i) normal mass hierarchy with the SMA or LMA scenario and a non-adiabatic 
H resonance, (ii) inverted mass hierarchy with the SMA or LMA scenario. In the 
antineutrino channels, significant earth matter effects are observed (sec. p.4|) only 
with (i) normal mass hierarchy with the LMA scenario, (ii) inverted mass hierarchy 
with LMA scenario and a non-adiabatic H resonance. 



7 Identification of the mass spectrum 

We show that with the 3z/ schemes which explain the solar and atmospheric neutrino 
data, one can make rather reliable predictions of the conversion effects for supernova 
neutrinos. The predictions differ for different schemes, which opens up the possibility 
of discriminating among them using the data from the neutrino bursts. 

We have studied the manifestation of the conversion effects in (i) The flavor com- 
position of the neutronization peak: the main effect here is the partial or complete 
change of the flavor of the peak. This can be established by the comparison of the 



44 





71 

P 


P 


I 


SMA 


normal 
inverted 


Pl 


1 




T A/TA 


IIUI llicLl 

inverted 


olli 1/0 


COS^ Oq 




V \J 


nui iiidi 

inverted 


kin^^^ po'^^a 1 

olli 5 ^^Vjo 


[sin^ 6^0, cos^ 6q] 


II 


SMA 


normal 
inverted 




1 

Ph 




T i\/rA 


noi mai 
inverted 


sin^ O^Ph 

olli (7q 


COS^ ^0 
IvUo UqIT u 




V \J 


noi mai 
inverted 


[\U,,\\cosHq] 

olli ^0 7 ^^vJo 1^0 


[sin^ 6^0, cos^ 6^0] 

[bin f70-r//,LOb UQjTffl 


III 


SMA 


normal 
inverted 


Pl 


1 
1 




LMA 


normal 
inverted 


sin^ ^0 
sin^ 6q 


COS^ ^0 

cos^ 6^0 




VO 


normal 
inverted 


[sin^ 610, cos^ 6q] 
[sin^ 6^0, cos^ ^0] 


[sin^ 6^0, cos^ 9q] 
[sin^ 00, cos^ ^0] 



Table 1: The values of the survival probabilities p and p for various scenarios. [x,y] 
indicates that the value of the survival probability lies between x and y. 

fluxes detected by charged current and neutral current interactions during the neu- 
tronization burst, (ii) Modifications of the z/g and Ue spectra: Here one expects an 
appearance of hard or composite and/or spectra due to mixing, instead of their 
original soft spectra. One can identify these effects of mixing by studying the average 
energies of spectra, the tails of spectra, and searching for the widening of the spectra, 
(iii) Earth matter effects: This can be done by detailed studies of the shapes of the 
energy spectra, or by the comparison of signals observed in different detectors. 

The final spectra of and Ue can be characterized by the values of the survival 
probabilities p and p respectively. The conversion probabilities depend significantly 
on the value of |?7e3p- We have divided the whole range of possible values of |?7e3p 
in three regions: I, II, and III, and have made definite predictions in each of these 
regions. Table |l] gives the values of the survival probabilities for the neutrino mass 
spectra under discussion. 

The qualitative features of the final neutrino fluxes for various neutrino mass 
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Tabic 2: The dependence of the final spectra on (1) the solution of the solar neutrino 
problem (2) the nature of hierarchy (normal or inverted) and (3) the value of |f/e3p. 
"Soft" in the z/g column refers to the original spectrum and "soft" in the column 
refers to the original Pe spectrum. "Hard" refers to the original spectrum. "Com- 
posite" implies that the final spectrum has both the soft and the hard components. 
In the Earth Matter Effects columns, a ^/ indicates the possibility of significant earth 
matter effects. 



46 



spectra, obtained in sec. are summarized in Table |^. Let us first make some 
general observations from the table. 

1. The complete disappearance of the Vg. neutronization peak (strong change of 
its flavor) and a pure hard spectrum of v^. during the cooling stage always appear 
together. Their observation will testify for l^/esl^ in the region I and normal mass 
hierarchy, irrespective of the solution of the solar neutrino problem. 

2. A pure hard z/g spectrum is the signature of the inverted mass hierarchy and 
|f/e3p in the region I, irrespective of the solar neutrino solution. 

3. A soft Ve spectrum can appear only in the schemes with the SMA solution of 
the solar neutrino problem. 

4. The observation of any earth matter effects will rule out the scenarios with the 
VO solution. 

5. The earth matter effects in z/g spectrum will be the signature of the LMA 
solution. 

6. If the earth matter effects are observed in the neutrino channels but not in 
the antineutrino channels, we either have the normal mass hierarchy with |?7e3p in 
regions II or III, or the inverted mass hierarchy. 

7. If the earth matter effects are significant in both the z/g and z/g channels, the 
solar neutrino solution is LMA and |t/e3p lies in regions II or III. 

Let us now systematically consider all the possible combinations of the observa- 
tions. In general, the "observations" consist of the following features of the spectra: 

• The flavor of the neutronization peak can be either mixed (z^g, ^x) or (almost) 
completely changed [y^. 

• The z/g spectrum at the cooling stage can be either hard or composite. 

• The z/g spectrum can be of all possible types: soft (unchanged), hard, composite. 

• The earth matter effects can be absent completely, or observed in one of channels 
(z/g or z/g), or observed in both channels. 

The following conclusions can be drawn from the combinations of the above ob- 
servations: 

(A) Completely changed flavor {ux) of the neutronization peak, hard z/g spectrum, 
and soft z/g spectrum. This configuration can be realized in only one case: normal 
hierarchy, SMA solution of the z/0-problem and |?7e3p in the region I. 

(B) Pure z/a; - neutronization peak, hard z/g spectrum, and composite z/g spectrum. 
This implies normal mass hierarchy, |f/g3p in the region I, and two possible Uq solu- 
tions, LMA or VO. The observation of earth matter effects can distinguish between 
these two possibilities. For the LMA solution, the earth effects will be present in the 
z/g channel and absent in the z/g channel. 

(C) Mixed (z/g, Ux) flavor neutronization peak, composite z/g spectrum, and soft 
z/g spectrum. This configuration can be realized in the three different cases: (i) SMA 
solution, normal mass hierarchy, and |t/e3p in the region II, (ii) SMA solution, normal 
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mass hierarchy, but |C/e3p in the region III, (iii) SMA solution, inverted hierarchy, 

|C/e3p in the region III. Clearly such a configTu^ation is the signature of the SMA 
solution of the i/0-problem. Earth matter effect is expected to be weaker in the case 



(D) Mixed (z/g, ^x) flavor neutronization peak, composite z/g spectrum, composite 
Pg spectrum. This combination of observations can be realized in many cases: (i) 
SMA solution, inverted hierarchy, |C/e3p in the region II, (ii-v) LMA solution, normal 
or inverted hierarchy, |f/e3p in the regions II or III, (vi-ix) VO solution, normal or 
inverted hierarchy, |t/e3p in the regions II or III. Studies of the earth effects can give 
an important criterion for discriminating between the above cases: they are absent 
in the VO cases (vi-ix) and should be observable in the i/g channel but not in Pg 
channel in the SMA case (i). In the LMA case, the earth effects arc expected both 
in the z/g and z/g channels. Here a further discrimination can be done between (ii-v) 
by studying the degree of compositness: in the case of normal hierarchy one expects 
stronger effects in the neutrino channel. 

(E) Mixed (i/g, flavor neutronization peak, composite fg spectrum, hard Pg 
spectrum. Such a combination of observations is realized with inverted hierarchy, 
|f/e3p in the region I, and SMA or LMA or VO solution. That is, the combination 
does not depend on the solution of the z/0-problem and turns out to be the signature 
of the inverted mass hierarchy. In the cases with SMA or LMA solutions the earth 
effect can be observed in the channel only. No earth effect is expected in the VO 
case. 

No other combination of the observables is realized in the mixing scenarios under 
disussion. 

Establishing the region in which |C/g3p lies is equivalant to giving bounds on the 
value of this mixing parameter, which correspond to the borders of the region. The 
borders depend on the precise density profile of the progenitor and can vary within a 
factor of 2 - 3. Thus, from SN data we can get bounds on |t/g3p with an uncertainty of 
a factor of 2 - 3. In principle the value of |^7g3p can also be measured (again within a 
factor of 2 -3) in the transition region where the effect (survival probability) depends 
substantially on Ues- 

The observations (A), (B), and (E) above establish |^7g3p in region I, thus giving a 
lower bound of |^7g3p > 3 • 10"''. The observations (C) and (D) establish |t/g3p to be 
in region II or III, thus giving an upper bound of |?7g3p < 3 • 10"''. The upper bound 
obtained in the case of the observations (C) and (D) is more than an order of mag- 
nitude stronger than the one the long baseline experiments are planning to achieve. 
The observation of the broadening of both the z/g and Pg spectra, or significant earth 
matter effects on both the z/g and Pg spectra are sufficient to establish this bound. 

A number of possible observations can rule out all the scenarios under discussion. 
The observation of (a) a pure i/g neutronization peak and (b) a soft z/g spectrum 
during the cooling stage is not possible with any of these neutrino mass spectra. The 
pairs of observations (a) pure z/,,- neutronization peak and a composite z/g spectrum, 
(b) mixed fiavour neutronization peak and a pure hard z/g spectrum, and (c) a hard 
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spectrum for both z/g and i>e are also fatal for all the scenarios. 

The final neutrino spectra can thus help in resolving three kinds of ambiguities 
that remain to be resolved with the current data: (i) the solution of the solar neutrino 
problem, (ii) the type of mass hierarchy (sign of Am^), and (iii) the value of |f/e3p- 

The implications of results of this paper will depend on when neutrino burst from 
the Galactic supernova will be detected. 

A number of results can be obtained from other experiments unrelated to the 
supernova neutrinos. Clearly, all the schemes considered above will be rejected if 
sterile neutrinos will be discovered in the oscillations of solar or atmospheric neutrinos, 
or if the LSND result will be confirmed. In this case, one will be forced to consider 
the schemes with sterile neutrinos 0. 

On the other hand, there is a good chance that within several years the existing 
and future experiments will allow us to identify one of solutions of the solar neutrino 
problem considered in this paper. This will significantly diminish the number of 
possible schemes and will allow us to further sharpen the predictions of the effects for 
supernova neutrinos. 

The purpose of this paper was to indicate the effects which can in principle testify 
for neutrino conversions and various features of the neutrino mass and flavor spectra. 
Clearly, further detailed studies are needed for specific detectors to clarify the de- 
tectability of the discussed effects and to conclude how far we can go in the program 
of identification of the neutrino mass spectrum. 
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